Topology of Smectic Order on Compact Substrates 
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Smectic orders on curved substrates can be described by differential forms of rank one (1-forrns), 
whose geometric meaning is the differential of the local phase field of density modulation. The 
exterior derivative of 1-form is the local dislocation density. Elastic deformations are described by 
superposition of exact differential forms. Applying this formalism to study smectic order on torus as 
well as on sphere, we find that both systems exhibit many topologically distinct low energy states, 
that can be characterized by two integer topological charges. The total number of low energy states 
scales as the square root of the substrate area. For smectic on a sphere, we also explore the motion 
of disclinations as possible low energy excitations, as well as its topological implications. 

PACS numbers: 61.30.-v,61.30.Jf,02.40.-k 



The problem of smectic order on curved substrate nat- 
urally arises in a variety of contexts, such as block copoly- 
mer films coated on a curved substrate [l], 0, d, H| , col- 
loidal particles immersed in a smectic liquid crystal @ 
with strong tangential boundary condition, smectic poly- 
mer vesicles, as well as Turing patterns on a sphere [6|. 
Last but not least, a flexible charged polymer adsorbed 
onto an oppositely charged curved surface Q may also 
form equal-distanced layer pattern on the surface. 

The main purpose of this letter is to study the low 
energy smectic states on curved substrates with mini- 
mal number of defects, which are experimentally most 
relevant at low temperature. We shall find that smectic 
order both on sphere and on torus exhibit many topologi- 
cally distinct, nearly degenerate, and minimally defected 
states. These states can be classified by two integer topo- 
logical charges, which are related to substrate topology 
as well as to disclination pattern. The total number of 
these states is proportional to VA/£o, where A is the sub- 
strate area, £q, the preferred layer spacing. Their energy 
differences scale sub-linearly with the system size. For 
smectic on sphere, by exchanging defects pairs in an ap- 
propriate manner, two integer charges can be systemat- 
ically changed. Interestingly, the corresponding process 
realizes the Euclidean algorithm for finding the greatest 
common divisor of two integers. Our results are comple- 
mentary to a recent work by Santangelo et al which 
addresses the energetic interplay between smectic order 
and substrate curvature. A more detailed version of our 
analysis is presented elsewhere Q . 

The local phases of density modulation of a smectic 
order at two nearby points x and x + dx are related by 
Q(x + Sx) « 9(x) + ip a (x)dx a . While the phase field 
is generically not globally well defined, its differential 
ip(x) = ip a dx a , called a I-form in modern differential 
geometry nevertheless still is. A simple geometric 
reasoning shows that the norm of ip is proportional to the 
reciprocal of layer spacing I: \ip\ — -J g al3 ip a i/j/3 = 2tt£~ 1 . 

The nonlinear strain W (x) = \ ( J^| 2 - l) measures 

the dilation/compression of smectic layers. This strain 
is identical to the one used in reference [l[ and reduces 
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FIG. 1: A): Two cycles on a torus. B) Smectic states on 
a torus with distinct global dislocation charges are shown as 
grid points in the upper-half (Ng,N$) plane. The curve (el- 
lipse) is the loci of vanishing strain w — 0, as determined by 
Eq. Q. States in the shaded region have low strain energy 
and are approximately degenerate. 



to the well known nonlinear strain for smectic order in 
flat space [ll], 03 • The total free energy of the system 
is the sum of the strain energy density Bw 2 /2 and the 
Frank free energy, appropriately generalized to curved 
space i[ '4 . 

The integral of ip along a loop 7 is the net phase change 
A0 as one walks around 7, i.e. the dislocation charges, 
multiplied by 2ir, enclosed by the loop. Using the Stoke's 
theorem this can in turn be related to the 2d integral over 
the region D enclosed by 7: 



A6 = 27TJV- 



dip, 



(1) 



where the two form dtp = (8 x ip y — 8 y ip x )dx A dy is the 
exterior differential of ip. dip therefore must be the dislo- 
cation density. Hence a dislocation-free smectic state is 
described by a 1-form satisfying dtp = 0. Such a 1-form 
is called closed. By contrast, a 1-form ip — d<3 that is the 
differential of a function O is called exact. 

An exact form is closed. If the phase field O is 
globally defined, we have ip = dO, and dip — d 2 Q — 
(8182 — 8281)0 dx 1 A dx 2 = 0, that is, the dislocation 
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density vanishes everywhere. This obvious result is called 
Poincare's lemma. Its converse holds in flat space: if 
dip = everywhere, the phase field is globally defined: 
ip = dQ, and the integral Eq. (fTJ) along arbitrary loop 
always vanishes. 

The converse of Poincare's lemma may fail in space 
with different topology [14[. De Rham's theorem of coho- 
mology [l5( | identifies the number of independent closed- 
but-not-exact forms with that of non-retractable loops 
on the manifold. A torus, for example, has two non- 
retractable loops, as shown in Fig. [TJ^. The most general 
closed 1-form on torus is given by 

il)N„N 4 =NedO + N+&t> + d$(6,<t>). (2) 

The integral of ip along loops 70 and 7^ give respectively 
2irN g and 2itN$. Eq. © is exact only if both Ng and 
N<f, vanish. In order for Eq. |3]| to describe a defects-free 
smectic state on torus, both Ng and must be integers. 
The exact form d$(6,<ft) in Eq. @ describes an elastic 
deformation of smectic layers relative to the state with 
d<& = 0. It can be ignored for a sufficiently thin torus. 

We define two states to be topologically equivalent if 
they can be brought into each other by continuous de- 
formation that do not break any smectic layer. The 
topological properties of a defects free smectic state on a 
curved substrate are therefore characterized by a closed- 
but-not-exact differential form. Similarly, any excita- 
tion/distortion of the smectic pattern that can not be re- 
laxed away by elastic deformation is a topological defects. 
According to this definition, therefore topological defects 
should be identified as the source of internal strain, a 
concept thoroughly explored in the continuous theory of 
defects in crystals [l6j]. We note, however, this definition 
of topological equivalence is more refined than the one 
used in the homotopy theory topological defects [13, EH • 
Consequently, the topological charges discussed in this 
work may not be truly topological invariants according 
to the homotopy theory. A more detailed discussion of 
this issue can be found in Q. 

Toroidal smectic states with different charges (Ng, N$) 
are clearly topologically distinct. Two integers (Ng, N^) 
share essential similarity with usual dislocation charges, 
but encode global properties of defects-free smectic states 
on a torus. Hence they shall be called global dislocation 
charges . Each defects-free smectic state can therefore 
be represented as a point in the half lattice of integers 
(N e ,N<f,) [2C|, as shown in Fig. 

The nonlinear strain associated with states Eq. (J5]) on 
a thin torus (with d$ set to zero) is given by: 

W - (2nRgf + (2*R+)* L W 

As shown in Fig. [TJ3, the equation w — traces out an 
ellipse in the (Ng,N$) plane. States satisfying < 
£ /^/A (shaded in Fig. QJ3) have the total strain energy 
bounded by Bi\j2, which is independent of the system 
size. A simple calculation Q also shows that the total 



Frank free energy scales sub-linearly with the system size. 
Therefore all these states are approximately degenerate 
for a large system. The total number of these low energy 
states is given by the area of the shaded region in Fig.[Tj3, 
and scales as the square root of the substrate area: 

Af(F B <l B f ) = 2^, (4) 

We now turn to smectic order on sphere. According 
to the Gauss-Bonnet-Poincare theorem, the total discli- 
nation charge of a smectic order on sphere must be two. 
Let us write ip = \ip\h, where h = n a dx a is the unit 1- 
form describing the smectic layer normal, while \tp\ is the 
norm of tp. Taking the exterior differential of ip, as well 
as using the Leibniz rule, we find dtp = (d\tp\) An + \ip\dh. 
Note that dh = (d\fi2 — d2ni)dx A dx 2 is a 2-form de- 
scribing bending deformation of the layer normal. Hence 
if there is no dislocation, dip = 0, and the layer spacing 
is constant, d\ip\ = 0, the bending deformation of the 
director field is strictly forbidden, dn = 0. This implies 
that the bending constant K3 is effectively infinity in a 
dislocation- free smectic (2lj . 

A spherical nematic in the limit of infinite bending 
rigidity is characterized by a one-parameter family of de- 
generate ground states, where four +1/2 disclinations sit 
on a great circle and form a rectangle with arbitrary as- 
pect ratio 22]. For a given bending-free nematic state, 
we can start from disclination cores and grow, layer by 
layer, a dislocation-free smectic pattern with equal layer 
spacing. We will, however, have to fine tune the layer 
spacing I so that the smectic pattern can be fit onto the 
sphere with given radius R. This fine tuning results in a 
small strain of order of £q/R, and a total strain energy 
Fb ~ B £q, which does not scale with the system size. 
Since the Frank free energy of all these states is the same 
by construction, the total free energy is approximately 
degenerate. 

Several numerical studies of spherical smectic orders 
formed by di-block copolymer films confined on sphere 
have been reported recently 0, 0|. By slowly an- 
nealing the system, three classes of low energy smec- 
tic states were found @: latitudinal states (23|, where 
all layers are circles of constant latitude, spiral states, 
where all smectic layers are spirals around the two poles, 
and quasi-baseball states, where four +1/2 disclinations 
are well separated on a great circle. All these states are 
dislocation-free. There may be, however, at most one 
layer termination at each disclination core. 

Consider a quasi-baseball state, shown in Fig.^V, with 
four disclinations sitting on the equator. We stereograph- 
ically project the pattern from the south pole to the com- 
plex plane, such that the equator is mapped into the unit 
circle. The image smectic pattern is shown in Fig. [5J3, 
with four defects located at (a, b, c, d) — (I, e 1 *, —I, e" 4 *) 
respectively. The segment ab of the unit circle intersects 
Ni/2 smectic layers. (A layer terminated at disclination 
a or b is counted as 1/2 intersection. ) Likewise, the seg- 
ment be of the unit circle intersects N2/2 smectic layers. 
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FIG. 2: (Color online) (A) A quasi-baseball state, with four 
+1/2 disclinations sitting on the equator. Each disclination 
core has one layer termination. Only two disclinations are 
visible here. (B) The same state stereographically projected 
onto the complex plane. The equator is mapped to the unit 
circle. The red contour encloses defects a, b, and intersects 
12 smectic layers in total. The blue contour encloses defects 
b, c and intersects 8 layers in total. Hence this quasi-baseball 
state is labeled by two integer charges (Ai = 12, A2 = 8). 

Note that by definition both Aq and A2 are non-negative. 
Without loss of generality we choose Nx < A2 . Two in- 
teger charges (Nx,N 2 ) completely determine the global 
properties of a quasi-baseball state. Latitudinal states 
are special cases of quasi-baseball states with Aq = [24| , 
while spiral states correspond to the regime Aq <C N 2 . 

It is clear from Fig. [2j\ that all smectic layers intersect 
the equator vertically. The arc lengths s a b and She are 
therefore given by 



s ab = R<t> = Nxt, 
from which we find 

N x _ (j> 

N 2 7T- < 



R(n 



nR 



Nx + N 2 



(5) 



where £ « £q is the layer spacing. The equilibrium value 
of the angle <p is therefore completely determined by the 
two charges Nx and N 2 - Large fluctuations of </> nec- 
essarily induces significant change of layer spacing, and 
therefore is energetically penalized. More importantly, 
Eq. ([5]) also shows that, in all low energy states, Aq + N 2 
is completely determined by the sphere radius R and the 
layer spacing I. Since Nx can take arbitrary integer value 
from to ttR/£, the total number of low energy smectic 
states on sphere is approximately given by 



M(F B <B£l) 




(6) 



i.e. also scales as the square root of the sphere area. 

To obtain an analytic description for quasi-baseball 
states, we invoke the Hodge decomposition theorem [l5[ 
of differential forms on compact manifold. The closed 
1-form ip can be expressed as the sum of the real part of 
an analytic form i\) c and an exact form da: 
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FIG. 3: (Color online) Starting from a state (Ni,N2) with 
Ai = 2, we twist the sphere counter-clockwise by n so as to 
exchange the defects pair a, b, and obtain a state (Ni,N2 + 
Ni). These two states are therefore topologically equivalent. 
This twist operation changes the layer spacing by a factor of 
Nx/N 2 . 



where ifj°(z) = f( z )dz, with f(z) — f(x + iy) a meromor- 
phic function, while a(x,y) is a smooth (but not ana- 
lytic) function on the z-plane. The topological properties 
of the smectic order on sphere is completely encoded in 
f(z)dz. The exact form da describes elastic deformation 
and should be chosen to minimize the total free energy. 

The complex meromorphic form t/j c that describes a 
spherical smectic with four disclinations at a, b, c, d can 
be shown to be M 



Ae 



y/ (z — a)(z — b)(z — c)(z — d) 



dz 



(8) 



As shown in Fig. OB, let us construct two loops 71 and 
72 that enclose the segments ab and be respectively. Inte- 
gration of ip along the loops 71 , 72 yield two topological 
charge Nx, N 2 that we defined above: 



(9) 



2TTN k =(b ip = Re(h i\) c , fc = l,2. 
J ik J ~tk 

Solving Eq. © for A and a we obtain 



Ae 



nN ? , 



ttNx 



k(oob£\ Af(sinf) 



J4>/2 



(10) 



4>(x, y) = Re ip c (z) + da(x, y), 



(7) 



where K(-) is complete elliptic integral of the first kind. 
Eqs. jS} and (fT0|) completely determine the meromorphic 
form tj) c for given charges Nx,N 2 . We note that the in- 
tegral of i/j c Eq. © defines the (multiple valued) inverse 
function of a doubly periodic elliptic function. Two in- 
teger charges Nx , N 2 are the real parts of two periods of 
this elliptic function [8j. 

Spiral state are especially interesting because of their 
closely bound defects pairs. As illustrated in Fig. [3l con- 
sider a spiral state with charges (Nx = 2,N 2 ), align the 
sphere so that one pair of defects a, b sit near the north 
pole. Let us fix the south pole and elastically twist the 
whole sphere, together with all smectic layers, around 
the z axis by an angle tt, so that the defects pair a, b at 
the north pole exchange their positions. The branch cut 
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(shown blue in Fig. [3j\) which connects a and d (invisible 
in Fig. [3]) before the twist operation, after the twist oper- 
ation, connects b and d instead, shown red in Fig.[3]B. The 
corresponding integer charge associated with this branch 
cut is N 2 before the twist, and becomes N 2 + Ni after 
the twist. Hence a counter-clockwise twist of the defects 
pair a, b by angle 7r leads to the following transformation 
of the topological charges 

(N 1 ,N 2 )-^(Nx 1 N 2 +Ni). (11) 

Obviously clockwise twist of the pair (a, b) leads to a 
different transformation: 

(N 1 ,N 2 )^(N U N 2 -N 1 ). (12) 

Therefore, if defects are allowed to move, states (iVi, N 2 ± 
Ni) are topologically identical to the state (N\, N 2 ). As 
long as A^i <C N 2 , this twist manipulation only induces 
small change in layer spacing and therefore costs little 
strain energy. Twist of closely bound defects pair is 
therefore a low energy excitation in spiral states. 

It is important to note that even if N\ and N 2 are com- 
parable, twist of defects pair is still topologically possible, 
but costs large strain energy. The two charges (Ni,N 2 ) 
of a quasi-baseball states are therefore stabilized by both 
topological and energetic barriers. Furthermore, defects 



twisting can be successively applied. Starting with an 
state with A^i < N 2 , we twist an appropriate defects pair 
so as to obtain new charges (iV{ = N\ , A^ — N 2 — N\ ) . 
If N[ < N2 still holds, we simply repeat the same 
twist operation. If, after the twist, we find N[ > N 2 , 
we twist a different pair so that we have a new state 
(N[ ,N2~N[). In every step, one of two integers (iVi , N 2 ) 
is reduced. This process therefore must end after finite 
steps, where one of the two integer charges vanishes. The 
final state is therefore a latitudinal state. More interest- 
ingly, this process of charge reduction is precisely the 
Euclidean algorithm for finding the greatest common di- 
visor gcd(Ni, N 2 ) between two integers (Ni,N 2 ). There- 
fore the nonvanishing charge of the final state must be 
gcd(Ni, N 2 ). Hence if disclinations are allowed to move, 
an arbitrary low energy smectic state on a sphere is topo- 
logically identical to certain latitudinal state, with only 
one nonzero integer topological charge. 
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